Electron states and electron Raman scattering in semiconductor double cylindrical quantum well wire
Munguía-Rodríguez M1, Betancourt-Riera Ri2, Betancourt-Riera Re2, †, , Riera R1, Nieto Jalil J M3
Departamento de Investigación en Física, Universidad de Sonora, Apartado Postal 5-88, 83190, Hermosillo, Sonora, México
Instituto Tecnológico de Hermosillo. Avenida Tecnológico S/N, Col. Sahuaro, 83170, Hermosillo, Sonora, México
Tecnológico de Monterrey-Campus Sonora Norte. Bulevar Enrique Mazón López No. 965, 83000, Hermosillo, Sonora, México

 

† Corresponding author. E-mail: rbriera@gmail.com

Abstract
Abstract

The differential cross section for an electron Raman scattering process in a semiconductor GaAs/AlGaAs double quantum well wire is calculated, and expressions for the electronic states are presented. The system is modeled by considering T = 0 K and also with a single parabolic conduction band, which is split into a subband system due to the confinement. The gain and differential cross-section for an electron Raman scattering process are obtained. In addition, the emission spectra for several scattering configurations are discussed, and interpretations of the singularities found in the spectra are given. The electron Raman scattering studied here can be used to provide direct information about the efficiency of the lasers.

1. Introduction

The Raman scattering is a useful technique to study the electronic structure of nanostructures. The analysis of the differential cross section of a Raman scattering process allows us to determine the subband structure of the nanostructured systems by a direct inspection of singularity positions in the spectra, considering the selection rules of transitions of the electron (or other quasiparticles) in the interaction with different polarizations of the incident and emitted light.[13] In general, the differential cross section of a Raman scattering process shows singularities related to the interband and intraband transitions, in our case, only the intraband transitions are considered.

The electron Raman scattering has been studied from a theoretical point of view in various systems. These studies considered the interband and intraband transitions with and without the participation of confinement phonons, also in the presence and absence of external electric and magnetic fields;[414] the bulk semiconductor, considering the presence of electric and magnetic fields, was also studied.[1517] Surface enhanced Raman scattering is one of the most powerful probing methods, which is non-destructive, real-time, and good at highly sensitive characterization.[18,19] Raman scattering studies of phonon modes are of fundamental importance for understanding the structural and optical properties, the surface states can lead to novel magnetic properties.[7,20,21] In the previous studies, the surface optical phonon mode, which is sensitive to the shell thickness, has been found in core/shell quantum dots and core–shell structured nanoparticles. Recently, it has been proposed to be realized in commonly used semiconductors, which suggests a promising approach to the development of semiconductor electronic devices.[21]

The development of several experimental techniques has allowed the nanometer-scale confinements of the band electrons in semiconductor materials to provide some varieties of quantum phenomena, such as low-dimensional electron states and modified dynamics of carriers in the systems. Thus, the study of electron Raman scattering remains as an interesting and fundamental issue.[22] Recently, the experimental research on the Raman scattering in the GaN semiconductor quantum well wire of cylindrical ring geometry has been reported.[2327] In particular, semiconductor double quantum well wire and multiple quantum wells have awakened great interest because these systems allow the development of light sources such as light emitting diodes, and laser diodes in a wide range of the electromagnetic spectrum from ultraviolet to terahertz. These systems were used because they offer the possibility to create an artificial three-level system[28,29] given the possibility of controlling the parameters of the quantum well wire such as width and height of the potential barrier. These devices are manufactured using different methods, and several types of materials, such as InGaAs/AlInAs and GaAs/AlGaAs.

In this work, we present a model of electron Raman scattering in double quantum well wire, where an electron is in the conduction band. The intraband Raman scattering processes can be qualitatively described in the following way: first, an electron in the conduction band absorbs a photon of incident radiation of energy, ħωl, resulting in an intersubband transition; then the electron emits a photon of secondary radiation of energy, ħωs, due to a new intersubband transition.

2. Model and electronic states

We must determine the stationary states of an electron in a double quantum well wire system of cylindrical geometry, where ra, rb, and rc are the radii and length L. The confinement potential (Vc) and the effective mass (μ) are given by

The envelope function approximation leads to the solution of Schrödinger equation

where

where u0 (r) is the electron Bloch function in the band; Jm and Nm are the Bessel functions of the first and the second kind of order m; Im and Km are the modified Bessel functions of the first and the second kind of order m. Here, we use the definitions of Bessel functions according to Abramowitz.[30] The electron energy is given by

where

where Er is the energy due to the spatial confinement. By considering the continuity of the function Ψ and the current density (1/μ)(∂Ψ/∂r) at the interface, we can calculate the constants A and B present in Eq. (2); also can determine energy levels Er, then

with

being xk = x (rk) and yk = y (rk), A1 can be obtained using the normalization condition, therefore,

where

In addition, Er is determined from the following transcendental equation:

where

As it can be observed, for each value of m, n energy values are obtained. The states are then described by the following quantum numbers: m = 0,1,2,…; n = 1,2,3,…; and kz (see Ref. [27]). We can say that the system is the combination of a quantum well wire of cylindrical disk geometry in the region where rra, which we call “the wire”, and a quantum well wire with a cylindrical ring geometry in the region where rbrrc, which we call “the shell”; while the barrier is the region where ra < r < rb.

3. Raman cross-section

The Raman differential cross-section in a volume per unit solid angle for incoming light of frequency ωl and scattered light of frequency ωs is given in Ref. [12], and it can be written as

where η(ωρ) is the refraction index as a function of the radiation with frequency ωρ and polarization unit vector is eρ. The subindex ρ = l(s) indicates the incident (secondary) radiation; c is the light velocity in a vacuum, V is the volume, and W(ωl,el,ωs,es) is the transition rate for the absorption and emission radiation calculated according to Fermi’s golden rule. In this case, two possibilities must be considered: (i) the electron absorbs the incident photon and it later emits the secondary radiation photon; or (ii) the electron emits the secondary radiation photon and later absorbs the incident photon.[28] Therefore,

where

in which |i⟩ and |f⟩ are the initial and final states with energies Ei and Ef, respectively; |a⟩ and |b⟩ are the intermediate states with energies Ea and Eb, respectively; and finally Γa and Γb are their respective lifetimes.[1]

The Hamiltonian operator for the radiation field Ĥρ, with frequency ωρ and polarization unit vector eρ, is commonly used in electron Raman scattering and resonant Raman scattering in nanostructures in the dipole approximation. This approximation is valid as long as the radiation wavelength is much larger than the radius of the wire. The effect of reduced dimensionality on the free carrier absorption and the selection rules for intersubband transition have been discussed in Ref. [28], where it is shown that the electron Raman scattering can only take place in the confinement direction, then

with μ the electron effective mass, and

Now, we must calculate the matrix elements appearing in Eq. (5). We use the condition k = 0 for the radiation field;[28] thus, by considering Eqs. (1) and (6), the following matrix elements can be obtained:

where and are the final and initial states of the transition, respectively, and

and μ0 is the electron free mass.

The explicit forms of Imm and IImm are very complicated and their expressions can be calculated following the method reported in Ref. [31]. As it can be observed, the emission or absorption of one photon can only take place between consecutive states of the quantum number m, according to the selection rules m″ = m′ ± 1, although there is no restriction for the quantum number n (see Eq. (7)).

In the initial state, we have an electron in a conduction subband and an incident radiation photon; and for the final state, we have an electron in the other conduction subband and a secondary radiation photon

For the intermediate states

Using Eqs. (7)–(9), we can obtain four different contributions to the differential cross section (see Ref. [27]), depending on the polarizations of the incident and emitted radiation, and finally we can rewrite Eq. (4) as

where

4. Discussion of the obtained results

In this case, we considered a cylindrical semiconductor double quantum wire grown in a GaAs/Al0.35Ga0.65 As matrix, and the following physical parameters[29] were used in this calculations: μ1 = 0.096μ0, μ2 = 0.0665μ0, and V0 = 0.30 eV, where μ0 is the free electron mass. Furthermore, in all cases, it is considered that Γa = Γb = Γf = 1 meV.

Figure 1 shows the electron state energy, Er of the double quantum well wire (solid lines), obtained using Eq. (3); where the values of ra and rc are fixed and rb takes values in the range of (ra, rc). In Figs. 1(a)1(c), ra = 6.0 nm and rc = 12.0 nm, ra = 6.0 nm and rc = 14.0 nm, ra = 8.0 nm and rc = 14.0 nm, respectively. The dashed lines correspond to the energies of a simple quantum well wire with radius ra. In this type of system, we can say that levels belong mainly to the wire or to the shell, but we must take into account that the quantum states actually belong to both, i.e., belong to the system; this is what we refer to as “coupling”.[32,33] Initially, in Fig. 1, we have a simple quantum well wire with radius rc, because rb = ra. As rb increases, the double quantum well wire is formed; such that the width of the shell, d = rcrb, diminishes and increases the width of the barrier, b = rbra. Finally, when rb = rc, we obtain a simple quantum well wire with radius ra. In Fig. 1(a), we can see that the electron state denoted by m = 0 and n = 1 is the ground state of the system, and is contained primarily in the wire, as it tends to a simple quantum well wire when the radius rb increases.[32] Besides, using the same argument we can determine that the electron state denoted by m = 0 and n = 2 is contained primarily in the shell, being the electron state of the shell with the lowest energy; but this is only true in the region rb > 6.6 nm; so this is the point where the system coupling begins to break. We can say that for 0 < rb < 6.6 nm, the system is in strong coupling. It can be observed that for rb ≈ 8.0 nm, the electron states corresponding to a shell are “taken out” from the system, due to the increase of the confinement of the electrons in this region. Moreover, we can see that the level denoted by m = 0 and n = 1 reaches the energy of a simple quantum well wire when rb ≈ 10.0 nm.

Fig. 1. Energy of electron states of the double quantum well wire (solid lines), where rb takes values in the range of (ra, rc), while (a) ra = 6.0 nm and rc = 12.0 nm, (b) ra = 6.0 nm and rc = 14.0 nm, and (c) ra = 8.0 nm and rc = 14.0 nm. The dashed lines correspond to the energies of a simple quantum well wire with radius ra.

If we compare Fig. 1(b) with Fig. 1(a), the rc value is increased, rc = 14.0 nm, being the immediate consequence of the increase in the number of electronic states “within” the system. Again, we have only one level of the wire, m = 0 and n = 1, which is logical as it remains the same value of ra = 6.0 nm, and again this is the lowest energy level. We can see that the electron state labeled as m = 0 and n = 2, cuts the level of the shell with the electron state, m = 1 and n = 1, approximately in rb = 8.0 nm, which determines the range where the coupling is strong. In this case, the strong coupling is in the range from b = 0 to b = 2.0 nm; that is larger than that obtained in Fig. 1(a). This implies that the increase of rc increases the size of the strong coupling region. Also, it is observed that the electronic states primarily associated with the shell, are organized in subbands determined by the quantum number n, see Fig. 1 of the Ref. [27]. In addition, we see that the electron state, which is denoted by m = 0 and n = 1, and associated mainly with the wire, reaches the energy of a simple quantum well wire when rb ≈ 10.0 nm. Comparing Fig. 1(a) with Fig. 1(b), we can conclude that increasing rc enlarges the strong coupling region, but practically does not change the point where the electron state denoted by m = 0 and n = 1 reaches the energy of a simple quantum well wire.

Moreover, in Fig. 1(c), the radii are set as ra = 8.0 nm and rc = 14.0 nm, which means that the radius of the wire is greater than that in Fig. 1(b) while the value of rc is the same. Comparing these two figures, we can find that the coupling increases because the barrier becomes smaller. Moreover, we can see that with the increase of the radius of the wire, the number of states that “fit” inside also increases, so there are primarily two energy levels associated with the wire. In these electron states mainly related to the wire, it is shown that the higher the energy is, the higher the coupling is, because the “barrier height” that must be overcome by the electrons is lower. When comparing these three figures in Fig. 1, we can see that the states with different values of m can be cut, while if it has the same value of m, it is not possible.

Figure 2 shows the energy of electron states, Er, for the double quantum well wire (black solid lines), where the barrier varies while the width of the shell is fixed, and the radii of the wire are ra = 6.0 nm and d = 6.0 nm, ra = 6.0 nm and d = 8.0 nm, and ra = 8.0 nm and d = 6.0 nm in Figs. 2(a)2(c), respectively. The dashed red lines correspond to the quantum state of a quantum well wire of cylindrical ring geometry with width d, while dashed blue lines correspond to the quantum state of a quantum well wire of cylindrical geometry with radius ra. In Figs. 2(a)2(c), we can see that all the electron states are corresponding with the band given by n = 1, as can be observed in Fig. 2 of Ref. [27]. Furthermore, in the three cases shown in Fig. 2, initially, the barrier does not exist, b = 0; this means that initially it has a simple quantum well wire with radius rc. By comparing Figs. 2(a)2(c), it can be deduced that the levels associated mainly to the wire do not cut each other, the same situation occurs with levels associated primarily to the shell, if they can cut the levels associated to different system components. In Fig. 2(a), it can be seen that as the barrier increases, the energy levels are adjusted to the energy levels corresponding to a shell (dashed red lines) and to a wire (dashed blue lines). For values of the barrier about 4.0 nm, we cannot find virtually any difference between the levels corresponding to the double quantum well wire and corresponding to the shell and to the wire, for any energy level. If we observe the levels associated with the shell, we can conclude that for higher energy, the decoupling is faster, because the “barrier height” that must be overcome by the electrons is lower. In the case, when electron states are associated primarily to the wire, the value of the barrier required for decoupling is approximately b = 1.5 nm, which is lower than that for the shell, so we can say that for the wire the decoupling is faster.

Fig. 2. Energy of electron states of the double quantum well wire (solid lines) where the barrier b = rbra varies, while ra and d = rcrb are maintaining fixed (a) ra = 6.0 nm and d = 6.0 nm, (b) ra = 6.0 nm and d = 8.0 nm, (c) ra = 8.0 nm and d = 6.0 nm. The dotted lines correspond to a quantum well wire of cylindrical ring geometry with width d, while dashed lines correspond to the energies of a simple quantum well wire with radius ra.

In Fig. 2(b), we increase the width of the shell at d = 8.0 nm, while maintaining the same internal radius, ra = 6.0 nm; this causes an increase in the number of levels mainly associated with the shell that are “within” the system, while it remains the number of associated states primarily to the wire. Comparing Fig. 2(b) with Fig. 2(a), we can see that the energy levels of double quantum wire are made virtually equal to the levels of the shell, from a width of about b = 2.5 nm. This means that when the width of the shell is increased, the electron states associated mainly to the shell are decoupled faster; so we can conclude that the coupling is weakened for the states of the shell. However, if we observe the blue dashed lines, it is noted that for the electron states associated primarily to the cable the decoupling requires a higher value of the barrier; in this case approximately b = 5.0 nm. Hence, we can conclude that in the case of states of the wire, the increase of d reinforces the coupling. Moreover, in Fig. 2(c) unlike in Fig. 2(a), we increase the value of ra to 8.0 nm, and the width of the shell, d = 6.0 nm remains constant. First, we note that two levels associated with the wire (dashed lines) appear, this is because if the radius is greater, more levels “within” the system can be found. Furthermore, it can be seen that the decoupling of the energy levels associated with the shell occurs to a value similar to the one seen in Fig. 2(a). On the other hand, in the case of the states associated to the wire, we can observe a coupling increase. Therefore, the coupling of the energy levels associated with the shell depends mostly on the width rather than the radius of the wire; while in the case of electron states associated with the wire, the coupling depends on the width and the radius of the shell and the wire, respectively.

Now, we have to analyze the differential cross section of the electron Raman scattering, which can explain the singularities presented in the emission and excitation spectra, respectively. In this case, we must note that equations (12) and (13) have resonant singularities, which are independent of the incident photon, or of the secondary photon when

where ρ = s is for the emission spectrum and ρ = l is for the excitation spectrum. While for Eqs. (10) and (11), we have the non-resonant singularities, which are dependent on the incident photon or the secondary photon for the emission spectrum and for the excitation spectrum, respectively, then

where a = 2 for and polarizations, and a = 0 for and polarizations. If the energy of the incident radiation is different to the resonance, this difference is named detuning.

Figure 3 shows the emission spectrum of a double quantum well wire with ra = 9.0 nm, rb = 12.0 nm, and rc = 17.0 nm, for incident photon with energy ħωl = 0.28 eV, where: (Fig. 3(a)), (Fig. 3(b)), (Fig. 3(c)). If we compare the contribution (Fig. 3(a)) with the contribution (Fig. 3(b)), we can see that the intensity of the resonant singularities is smaller, moreover it can also be noted that the position of singularities change. This is because the two figures correspond to two different polarizations of secondary radiation, , which affects the selection rules of m ± 1 = m′. As we can see, the position of singularities in the graph is related to the “gap” between the two energy levels involved, while the peak intensity is mainly related to the probability of the transition and the detuning (see Ref. [29]), if the incident radiation (ħωl) is the same. When comparing Fig. 2(c) with Figs. 2(a) and 2(b), it is easy to see that the intensity of the resonant singularity is higher by almost an order of magnitude, that being the reason for the main contribution to the differential cross section. Also, we can see that the resonant singularities in polarizations and are the same, this is because the secondary radiation fulfills in both cases the identical selection rule.

Fig. 3. The emission spectrum for an electron Raman scattering process in double quantum well wire with ra = 9.0 nm, rb = 12.0 nm, and rc = 17.0 nm; for incident photon with energy ħωl = 0.28 eV and the polarization: (a) ; (b) ; (c) .

Moreover, figure 4 shows the emission spectrum for an electron Raman scattering process in double quantum well wire for incident photon energy ħωl = 0.28 eV and polarization . In Fig. 4(a), the dashed line corresponds to ra = 8.0 nm, rb = 10.0 nm, and rc = 15.0 nm, while the solid line corresponds to ra = 9.0 nm, rb = 10.0 nm, and rc = 15.0 nm. In Fig. 4(a), the dashed line corresponds to ra = 8.0 nm, rb = 10.0 nm, and rc = 16.0 nm, while the solid line corresponds to ra = 8.0 nm, rb = 12.0 nm, and rc = 17.0 nm. In Fig. 4(b), the dashed line has been multiplied by 10. As it can be seen in Fig. 4(a) that the increase of the radius of the wire, ra, causes an increase in the number of singularities and position changing of singularities, which we can already see in the spectrum, and these singularities move to a lower energy position. This is because the increase of ra increases the number of quantum states “within” the system, especially those mainly related to the wire, while the gap between different subbands decreases (see Eq. (14)). This is corroborated if we observe the resonant singularity denoted by ħωs (1,1,0,1) in both figures. In addition, in Fig. 4(b), we maintain a constant radius of the wire, ra = 8.0 nm, and change the width of the shell and the barrier. However, the position of the singularity ħωs(1,1,0,1), remains constant (also, note the dashed line in Fig. 4(a); this confirms that the quantum states m = 1, n = 1 and m = 1, n = 1 belong mainly to the wire. If we look at the singularity ħωs (1,2,0,1), this moves to a position of greater energy, which is because the quantum state labeled as m = 1, n = 2, is mainly related to the shell. Furthermore, we can compare the dashed line in Fig. 4(a) with the dashed line in Fig. 4(b). In this case, the width of the shell is increased by rc, increasing and it causes a similar effect to the one described in Fig. 4(a), but mainly affecting quantum states of the shell, allowing the occurrence of singularities related to the ones with the energy levels of the shell. Whereby, we clarify that the increase of ra mainly affects the energy levels belonging to the wire, while the increase of rc primarily affects the energy levels belonging to the shell. Finally, the dashed line in Fig. 4(a) and the solid line in Fig. 4(b) can be compared. We can see that the increment of the barrier, b, keeping a constant width, d, only causes a small change in the position of the resonant peak, which corresponds with that observed in Fig. 2, but affects the intensity of the singularities.

Fig. 4. The emission spectrum for an electron Raman scattering process in double quantum well wire for incident photon energy ħωl = 0.28 eV and polarization , where: (a) the dashed line corresponds to ra = 8.0 nm, rb = 10.0 nm, and rc = 15.0 nm, while the solid line corresponds to ra = 9.0 nm, rb = 10.0 nm, and rc = 15.0 nm; (b) the dashed line corresponds to ra = 8.0 nm, rb = 10.0 nm, and rc = 16.0 nm, while the solid line corresponds to ra = 8.0 nm, rb = 12.0 nm, and rc = 17.0 nm. In panel (b), the dashed line has multiplied by 10.

The change in the intensity of differential cross section for an electron Raman scattering process in semiconductor double quantum well wire is related to the coupling between the wire and the shell. This is evidently shown in Fig. 5, in which a transition rate that describes the following process is shown: first, an electron in the conduction band absorbs a photon of incident radiation, resulting in a transition from initial state m″ = 0, n″ = 1 to intermediate state m′ = 1, n′ = 1; then, the electron emits a photon of secondary radiation due to a new transition, from intermediate state m′ = 1, n′ = 1 to final state m = 0, n = 2. It is important to note that the net Raman gain in this case is directly proportional to the transition rate, this is |T (0,1)T+ (1,0)|2. Then, figure 5 shows the rate transition to a wire for a double quantum well, where the barrier b varies while the fixed values of ra and d remain. In addition, we must note that the solid line corresponds to a double quantum well where ra = 8.0 nm and d = 6.0 nm, being the dash line the one that corresponds to a double quantum well where ra = 8.0 nm and d = 8.0 nm. These curves are much like those obtained from an asymmetric multiple quantum well[32] and different from those obtained for a step-quantum well;[29] this is because the double quantum well wire can be decoupled as the asymmetric multiple quantum well, but in the present case the coupling is weakened much faster. Summarizing, the coupling in a double quantum well wire with cylindrical symmetry is weaker than an asymmetrical double quantum well. In addition, it is noted that any change in the system parameters causes strong changes in the transition rate, making this system very difficult to optimize. However, although the coupling is weak in this system, the transition rate is larger than the one obtained for the asymmetrical double quantum well, observed in Fig. 6 in Ref. [32], which is due to the symmetry of the system. We should clarify that when we refer to the coupling, we are referring to the interaction between the parts of the active region of the system, which in this case are the wire and the shell. This interaction affects the electronic states of the system, so that the electronic states are a mix between wire electron states and shell electron states.

Fig. 5. The transition rate, |T (0,1)T+ (1,0)|2, of a double quantum well wire, where the barrier b = rbra varies while ra and d = rcrb are fixed. The solid line corresponds to ra = 8.0 nm and d = 6.0 nm, while the dash line corresponds to ra = 8.0 nm and d = 8.0 nm.
5. Conclusions

In this work, we presented the formulas of calculating the electron Raman scattering intensities for the semiconductor double quantum well wire with a cylindrical cross-section, grown in a GaAs/Al0.35Ga0.65 As matrix. Moreover, we demonstrated that the structure of differential cross-section as given in figures, provides a transparent understanding of the energy subband structure. Furthermore, we showed that any change in the system configuration leads to changes in the intensity and the frequency of the secondary radiation, as well as the net Raman gain. This effectively shows that Raman spectroscopy is a powerful probe for semiconductor nanostructures due to the material-selective resonant enhancement. In addition, we can see that in this work we have applied a simplified model for the electron structure of the system. In a more realistic case, we should consider a coupled band structure by using models like the Luttinger–Kohn or the Kane model. Anyway, this band model allows us to analyze the main features of this system easily.

Reference
1Betancourt-Riera RRiera RMarín J LRosas R A2004Electron Raman Scattering in Nanostructures, Encyclopedia of Nanoscience and NanotechnologyStevenson RanchAmerican Scientific PublishersVol.3101
2Kushwaha M S 2001 Surf. Sci. Rep. 41 1
3Riera RMarín J LRosas R A2001Optical Properties and Impurity States in Nanostructured Materials, Handbook of Advanced Electronic and Photonic DevicesNew YorkAcademicVol. 6Chapter 6
4Betancourt-Riera RBetancourt-Riera RNieto J J MRiera R 2013 Physica 410 126
5Betancourt-Riera RBetancourt-Riera RRosas RRiera R 2012 J. Comput. Theor. Nanosci. 9 1
6Lu FLiu C HGuo Z L 2012 Physica 407 165
7Zhai L XWang YLiu J J 2011 J. Appl. Phys. 110 043701
8Zhong Q HYi X H 2010 Superlatt. Microstruct. 47 723
9Zhong Q HLiu C HZhang Y QSun H C 2008 Phys. Lett. 372 2103
10Scheinert MSigg HTsujino S 2007 Appl. Phys. Lett. 91 131108
11Ismailov T GMehdiyev B H 2006 Physica 31 72
12Bergues J MBetancourt-Riera RRiera RMarin J L 2000 J. Phys.: Condens. Matter 12 7983
13Bergues J MBetancourt-Riera RMarin J LRiera R1996Phys. Low-Dimens. Struct.7/881
14Riera RComas FTrallero-Giner CPavlov S T 1988 Phys. Status Solid 148 533
15Comas FTrallero G CPerez-Alvarez R 1986 J. Phys. C: Solid State Phys. 19 6479
16Trallero-Giner CRuf TCardona M 1990 Phys. Rev. 41 3028
17González de la Cruz GTrallero-Giner C 1998 Phys. Rev. 58 9104
18Wang M LZhang C XWu Z LJing X LXu H J 2014 Chin. Phys. 23 067802
19Jiang S MWu D JWu X WLiu X J 2014 Chin. Phys. 23 047807
20Wang MTian YZhang J MGuo C FZhang X ZLiu Q 2014 Chin. Phys. 23 087803
21Shi L KLou W K 2014 Chin. Phys. Lett. 31 067304
22Zhong Q HSun Y T 2011 Thin Solid Films 519 8178
23Wang JDemangeot FPechou RPonchet ACros ADaudin B 2012 Phys. Rev. 85 155432
24Livneh TZhang J PCheng G SMoskovits M 2006 Phys. Rev. 74 035320
25Laneuville VDemangeot FPchou RSalles PPonchet AJacopin GRigutti Lde Luna Bugallo ATchernycheva MJulien F HMarch KZagonel L FSongmuang R 2011 Phys. Rev. 83 115417
26Liu H LChen C CChia CYeh C CChen C HYu M YKeller SDenBaars S P 2001 Chem. Phys. Lett. 345 245
27Betancourt-Riera RBetancourt-Riera RNieto J J MRiera R 2015 Chin. Phys. 24 117302
28Betancourt-Riera RNieto J J MRiera RBetancourt-Riera RRosas R 2008 J. Phys.: Condens. Matter 20 045203
29Ferrer-Moreno L ABetancourt-Riera RBetancourt-Riera RRiera R2015Physica B47787
30Abramowitz MStegun I A1972Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical TablesNew YorkDover Publications10.1016/j.physb.2015.08.024
31Bergues J MRiera RComas FTrallero-Giner C 1995 J. Phys.: Condens. Matter 7 7273
32Betancourt-Riera RRosas RMarín-Enriquez IRiera RMarin J L 2005 J. Phys.: Condens. Matter 17 4451
33Betancourt-Riera RRiera RRosas RNieto J M 2008 Physica 40 785